Modular curves, C*-algebras, and chaotic 
cosmology 
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We make some brief remarks on the relation of the mixmaster universe model 
of chaotic cosmology to the geometry of modular curves and to noncommu- 
tative geometry. We show that the full dynamics of the mixmaster universe 
is equivalent to the geodesic flow on the modular curve Xp t 2 )- We then con- 
sider a special class of solutions, with bounded number of cycles in each Kas- 
ner era, and describe their dynamical properties (invariant density, Lyapunov 
exponent, topological pressure). We relate these properties to the noncom- 
mutative geometry of a moduli space of such solutions, which is given by a 
Cuntz-Krieger C*-algebra. 

1 Modular curves 

Let G be a finite index subgroup of r = PGL(2,Z), and let Xq denote 
the quotient Xq = G\M 2 , where H 2 is the 2-dimensional real hyperbolic 
plane, namely the upper half plane {z G C : 5z > 0} with the metric ds 2 — 
\dz\ 2 / (^sz) 2 . Equivalently, we identify H 2 with the Poincare disk {z : \z\ < 1} 
with the metric ds 2 = 4|dz| 2 /(l — |z| 2 ) 2 . 

Let P denote the coset space P = r/G. We can write the quotient Xq 
equivalently as X G = T\(H 2 x P). The quotient space Xq has the structure 
of a non-compact Riemann surface, which can be compactified by adding cusp 
points at infinity: 



In particular, we consider the congruence subgroups G = i~b(p), with p a 
prime, given by matrices 



satisfying c = mod p. In fact, for our purposes, we are especially interested 
in the case p — 2. 




(1) 
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1.1 Shift operator and dynamics 

If we consider the boundary P 1 (R) of H 2 , the arguments given in [5] and [TU] 
show that the quotient I n \(P 1 (M) x P), better interpreted as a "noncommuta- 
tive space" , gives rise to a compactification of the modular curve Xq with a 
structure richer than the ordinary algebro-geometric compactification by cusp 
points. 

In 9 this was described in terms of the following dynamical system, gen- 
eralizing the classical Gauss shift of the continued fraction expansion: 

T : [0, 1] x P -> [0, 1] x P 



T(x,t) 



1 0) 



(2) 



In fact, the quotient space of P 1 (M) x P by the PGL(2,Z) can be identified 
with the space of orbits of the dynamical system T on [0, 1] x P. 

We use the notation x — [ki, k%, . . . , k n , . . .] for the continued fraction ex- 
pansion of the point x £ [0, 1], and we denote by p n {x) / q n {x) the convergents 
of the continued fraction, with p n (x) and q n (x) the successive numerators and 
denominators. It is then easy to verify that the shift T acts on x by shifting 
the continued fraction expansion, T[k±, fe, ■ ■ ■ , k n , . . .] = [fo, Aft, . . . , k n +i, . . .]. 
The element g n (x)~ 1 , with 

9n(x) = ( Pk -f\ p f\)er, 

yny ' \q k -l(x) q k (x) J 

acts on [0, 1] x P as T n . 
The Lyapunov exponent 

X(x) := lim -log\(T n )'(x)\ 

n — >oc Ji 

= 2 lim — log q n (x) 

n — >oo ji 

measures the exponential rate of divergence of nearby orbits, hence it provides 
a measure of how chaotic the dynamics is. 

For the shift of the continued fraction expansion, the results of ^3] show 
that A(ir) = 7r 2 /(61og2) = A almost everywhere with respect to the Lebesgue 
measure on [0, 1] and counting measure on P. On the other hand, the Lyapunov 
exponent takes all values X(x) € [Ao,oo). The unit interval correspondingly 
splits as a union of T-invariant level sets of A (Lyapunov spectrum) of varying 
Hausdorff dimension, plus an exceptional set where the limit defining A does 
not converge. 
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2 Mixmaster universe 

An important problem in cosmology is understanding how anisotropy in the 
early universe affects the long time evolution of space-time. This problem is 
relevant to the study of the beginning of galaxy formation and in relating 
the anisotropy of the background radiation to the appearance of the universe 
today. 

We follow [5] (cf. also ^21) for a brief summary of anisotropic and chaotic 
cosmology. The simplest significant cosmological model that presents strong 
anisotropic properties is given by the Kasner metric 

ds 2 = -dt 2 + t 2pi dx 2 + t 2p2 dy 2 + t 2p3 dz 2 , (3) 

where the exponents ft are constants satisfying J^P' = 1 = J2iPi- Notice 
that, for pi = dloggu/dlogg, the first constraint J2%Pi = 1 is just the condi- 
tion that log gij — 2aSij + for a traceless f3, while the second constraint 
— 1 amounts to the condition that, in the Einstein equations written 
in terms of a and , 

= — (t oq I 1 
dt J 3 y 16tt V dt 

p ~s a d ( e3a dPij\ ( 1 

6 Jt\ e -dT) = ^\ r ^--^ Tkk 

the term T 00 is negligible with respect to the term (dfiij /dt) 2 /16tt, which 
is the "effective energy density" of the anisotropic motion of empty space, 
contributing together with a matter term to the Hubble constant. 

Around 1970, Belinsky, Khalatnikov, and Lifshitz introduced a cosmolog- 
ical model (mixmaster universe) where they allowed the exponents pi of the 
Kasner metric to depend on a parameter u, 

Pi 




1 + U + U 2 

P2 = (4) 

_ m(1+m) 
P3 ~ 1 + u+u 2 

Since for fixed u the model is given by a Kasner space-time, the behavior 
of this universe can be approximated for certain large intervals of time by a 
Kasner metric. In fact, the evolution is divided into Kasner eras and each era 
into cycles. During each era the mixmaster universe goes through a volume 
compression. Instead of resulting in a collapse, as with the Kasner metric, high 
negative curvature develops resulting in a bounce (transition to a new era) 
which starts again a behavior approximated by a Kasner metric, but with 
a different value of the parameter u. Within each era, most of the volume 
compression is due to the scale factors along one of the space axes, while the 
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other scale factors alternate between phases of contraction and expansion. 
These alternating phases separate cycles within each era. 
More precisely, we are considering a metric 

ds 2 = -dt 2 + a{t)dx 2 + b(t)dy 2 + c(t)dz 2 , (5) 

generalizing the Kasner metric ©■ We still require that © admits SO(3) 
symmetry on the space-like hypersurfaces, and that it presents a singularity 
at t — -> 0. In terms of logarithmic time dQ = the mixmaster universe 

& abc ' 

model of Belinsky, Khalatnikov, and Lifshitz admits a discretization with the 
following properties: 

1. The time evolution is divided in Kasner eras [!?„, i7 n+ i], for n £ Z. At the 
beginning of each era we have a corresponding discrete value of the parameter 
u n > 1 in (@J. 

2. Each era, where the parameter u decreases with growing f2, can be sub- 
divided in cycles corresponding to the discrete steps u n , u n — 1, u n — 2, etc. 
A change u — > u — 1 corresponds, after acting with the permutation (12) (3) 
on the space coordinates, to changing u to —it, hence replacing contraction 
with expansion and conversely. Within each cycle the space-time metric is 
approximated by the Kasner metric J21 with the exponents pi in Q with a 
fixed value of u = u n — k > 1 . 

3. An era ends when, after a number of cycles, the parameter u n falls in the 
range < u n < 1. Then the bouncing is given by the transition u — > 1/u which 
starts a new series of cycles with new Kasner parameters and a permutation 
(1)(23) of the space axis, in order to have again pi < p 2 < P3 as in 

Thus, the transition formula relating the values u n and u n+ \ of two suc- 
cessive Kasner eras is 

1 



which is exactly the shift of the continued fraction expansion, Tx = 1/x— [1/x] , 
with x n+ i = Tx n and u n — l/x n . 



3 Geodesies and universes 

The previous observation is the key to a geometric description of solutions 
of the mixmaster universe in terms of geodesies on a modular curve (Manin- 
Marcolli 0): 

Theorem 3.1 Consider the modular curve Xr (2)- Each infinite geodesic 7 
on Xp (2) not ending at cusps determines a mixmaster universe. 

Proof. An infinite geodesic on Xp ^ 2 ) is the image under the quotient map 
Ti r : H 2 x P -> T\(H 2 x P) X G , 
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where r = PGL(2,Z), G = T (2), and P = FjG ^ P^Fa) = {0,1, oo}, 
of an infinite geodesic on H 2 x P with ends on P 1 (R) x P. We consider the 
elements of P 1 (F 2 ) as labels assigned to the three space axes, according to the 
identification 

= [0 : 1] i-> z 

00 = [1 : 0] i-> y (6) 

1 = [1 : 1] ^ x. 

Any geodesic not ending at cups can be coded in terms of data (u>~ , uj + , s) , 
where (lu^, s) are the endpoints in P 1 (K) x {s}, s G P, with u~ G (— oo, —1] 
and lj + G [0, 1]. In terms of the continued fraction expansion, we can write 

LO + = [ko, kl,...kr, kr+l, ■ ■ ■] 



The shift acts on these data by 



T(u-,a) = 



1 

1 



w- + [l/(j+]'\ 1 y 

Geodesies on Xp ^ 2 ) can be identified with the orbits of T on the set of data 
(w,s). 

The data (w, s) determine a mixmaster universe, with the k n = [w„] = 
[l/x n ] in the Kasner eras, and with the transition between subsequent Kasner 
eras given by x n+ \ — Tx n G [0, 1] and by the permutation of axes induced by 
the transformation 

-k n 1 
1 / 

acting on P 1 ^). It is easy to verify that, in fact, this acts as the permutation 
i — ► oo, 1 i — ► 1, oo i — ► Q, if k n is even, and i— > oo, 1 i— > 0, oo i— > 1 if k n 
is odd, that is, after the identification JJjJl, as the permutation (1)(23) of the 
space axes (x, y, z), if k n is even, or as the product of the permutations (12)(3) 
and (1)(23) if k n is odd. This is precisely what is obtained in the mixmaster 
universe model by the repeated series of cycles within a Kasner era followed 
by the transition to the next era. 

Data (to,s) and T m (Lj,s), m G Z, determine the same solution up to a 
different choice of the initial time. 

There is an additional time-symmetry in this model of the evolution of 
mixmaster universes (c/. In fact, there is an additional parameter S n in 
the system, which measures the initial amplitude of each cycle. It is shown in 
that this is governed by the evolution of a parameter 

5 n+ i(l + u n ) 
Vn = ~ i — x 

1 - On+l 
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which is subject to the transformation across cycles v n+ i — [u n ] + v n . By 
setting y n — v^ 1 we obtain 

1 



Vn+l 



(yn + [1/Xn]) ' 



hence we see that we can interpret the evolution determined by the data 
(uj + ,uj~,s) with the shift T either as giving the complete evolution of the 
u-parameter towards and away from the cosmological singularity, or as giving 
the simultaneous evolution of the two parameters (u, v) while approaching the 
cosmological singularity. 

This in turn determines the complete evolution of the parameters (it, 5, i?), 
where fl n is the starting time of each cycle. For the explicit recursion J7 n +i = 
O n+ i(f2 n ,x n ,y n ) see 

Notice that, unlike the description of the full mixmaster dynamics given, 
for instance, in J2J , we include the alternation of the space axes at the end of 
cycles and eras as part of the dynamics, which is precisely what determines the 
choice of the congruence subgroup. This also introduces a slight modification 
of some of the invariants associated to the mixmaster dynamics. For instance, 
it is proved in [H] that there is a unique T-invariant measure on [0, 1] x P, 
which is given by the Gauss density on [0, 1] and the counting measure 5 on 
P: 

, / s S(s)dx . . 

rf/i(2; ' s) = 31og(2)(l + x)' (7) 
which reduces to the Gauss density for the shift of the continued fraction on 
[0, 1] when integrated in the P direction. In particular, as observed in the 
form {JJ of the invariant measure implies that the alternation of the space 
axes is uniform over the time evolution, namely the three axes provide the 
scale factor responsible for volume compression with equal frequencies. 



4 Controlled pulse universes 

The interpretation of solutions in terms of geodesies provides a natural way 
to single out and study certain special classes of solutions on the basis of 
their geometric properties. Physically, such special classes of solutions exhibit 
different behaviors approaching the cosmological singularity. 

For instance, the data (lo + ,s) corresponding to an eventually periodic 
sequence koki . . . k m . . . of some period ciq . . . ai correspond to those geodesies 
on Xp / 2 ) that asymptotically wind around the closed geodesic identified with 
the doubly infinite sequence . . . eto . . . aeao . . . ag . . .. Physically, these universes 
exhibit a pattern of cycles that recurs periodically after a finite number of 
Kasner eras. 

In the following we concentrate on another special class of solutions, which 
we call controlled pulse universes. These are the mixmaster universes for which 
there is a fixed upper bound N to the number of cycles in each Kasner era. 
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In terms of the continued fraction description, these solutions correspond 
to data s) with ui + in the Hensley Cantor set E^ C [0, 1]. The set is 
given by all points in [0, 1] with all digits in the continued fraction expansion 
bounded by N {cf. 6 ). In more geometric terms, one considers geodesies on 
the modular curve Xr (2) that wander only a finite distance into the cusps. 



4.1 Dynamical properties 

It is well known that a very effective technique for the study of dynamical 
properties such as topological pressure and invariant densities is given by 
transfer operator methods. These have already been applied successfully to 
the case of the mixmaster universe, cf. In our setting, for the full mix- 
master dynamics that includes alternation of space axes, the Perron- Frobcnius 
operator for the shift J3J) is of the form 



(x + k)! 3 \x + k'' \1 k 



This yields the density of the invariant measure 10 satisfying £2/ = /• 
The top eigenvalue rip of Cp is related to the topological pressure by 77^ = 
exp(P(/3)). This can be estimated numerically, using the technique of pQ and 
the integral kernel operator representation of §1.3 of 0. 

We now restrict our attention to the case of controlled pulse universes. 
Since these form sets of measure zero in the measure (JJJ , they support excep- 
tional values of such dynamical invariants as Lyapunov exponent, topological 
pressure, entropy. In fact, for a fixed bound N on the number of cycles per 
era, we are considering the dynamical system (|2} 

T : E N xP^ E N x P. 

For this map, the Perron-Frobenius operator is of the form 

(w)^ S )=e^/(^(;j)-). (8) 

It is proved in [10) that this operator still has a unique invariant measure ynjv, 
whose density satisfies £ 2 dim H (E N ),Nf = f, with 

the Hausdorff dimension of the Cantor set En- Moreover, the top eigenvalue 
r]p of C/3 t N is related to the Lyapunov exponent by 

d 



for fi N -almost all x £ E N , cf. [TU] . 
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5 Non-commutative spaces 

A consequence of this characterization of the time evolution in terms of the 
dynamical system @ is that we can study global properties of suitable moduli 
spaces of mixmaster universes. For instance, the moduli space for time evolu- 
tions of the u-parameter approaching the cosmological singularity as fl — > oo 
is given by the quotient of [0, 1] x P by the action of the shift T. Similarly, 
the moduli spaces that correspond to controlled pulse universes are the quo- 
tients of En x P by the action of the shift T. It is easy to see that such 
quotients are not well behaved as a topological spaces, which makes it dif- 
ficult to study their global properties in the context of classical geometry. 
However, non-commutative geometry in the sense of Connes 0] provides the 
correct framework for the study of such spaces. The occurrence in physics of 
non-commutative spaces as moduli spaces is not a new phenomenon. A well 
known example is the moduli space of Penrose tilings (cf. 0]), which plays an 
important role in the mathematical theory of quasi-crystals. 

We consider here only the case of controlled pulse universes. In this case, 
the dynamical system T is a subshift of finite type which can be described by 
the Markov partition 

A N = {((fc, t), (£, s))\U k , t C T(Ui, s )}, 

for k, £ G {1, . . . , N} : and s, t G P, with sets U k ,t = Uk x {*}, where U k C E N 
are the clopen subsets where the local inverses of T are defined, 



1 



1 



1' k 



DE N . 



This Markov partition determines a matrix An, with entries (A/v)/ct,fe = 1 if 
Uk,t C T(Ue tS ) and zero otherwise. 

Lemma 5.1 The 3x3 submatrices A k g = (^4(fe ) t) ) (^ lS )) S) tgp of the matrix An 
are of the form 



A u = I 





(001\ 


M x = 


10 


< 


vioo; 




(001\ 


M 2 = 


10 




voioy 



2 m 



The matrix An is irreducible and aperiodic. 

Proof. The condition U kl t C T(£/g s ) is equivalently written as the condition 
that 

'0 r 
1 £ 



t. 
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We then proceed as in the proof of Theorem 13 . II and notice that the transfor- 
mation 

'or 
i t t 

acts on P 1 (F2), under the identification ©, as the permutation i— » oo, 1 i— * 1, 
oo I— > 0, when £ is even, and oo i— * 0, i— » 1, 1 1— ► oo if £ is odd. 

Irreducibility of Am means that the corresponding directed graph is 
strongly connected, namely any two vertices are connected by an oriented 
path of edges. Since the matrix An has the form 



.4 



/ Mi Mi Mi 
M 2 M 2 M 2 
Mi Mi Mi 

V : 



Mi\ 

M 2 

Mi 

: / 



irreducibility follows from the irreducibility of Ai , which corresponds to the 
directed graph illustrated in the Figure. This graph also shows that the matrix 



Am is aperiodic. In fact, the period is defined as the gcd of the lengths of the 
closed directed paths and the matrix is called aperiodic when the period is 
equal to one. 

As a non-commutative space associated to the Markov partition we con- 
sider the Cuntz-Krieger C*-algebra Oa n (cf. 0), which is the universal C*- 
algebra generated by partial isometries Skt satisfying the relations 

22 SmSm = i, 
(fc,*) 

(k,t) 
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It is a well known fact that the structure of this C*-algebra reflect proper- 
ties of the dynamics of the shift T. For example, one can recover the Bowen- 
Franks invariant of the dynamical system T from the X-theory of a Cuntz- 
Krieger C*-algebra Oa (cf. 0). Moreover, in our set of examples, information 
on the dynamical properties of the shift T and the Perron-Frobenius operator 
JHJ can be derived from the KMS states for the C* -algebras Oa n with respect 
to suitable time evolutions. 

5.1 Time evolution and KMS states 

Recall that a state tp on a unital C*-algebra A is a continuous linear functional 
f : A — > C satisfying ip(a*a) > and <p(l) = 1. Let at be an action of 
M on A by automorphisms. A state ip satisfies the KMS condition at inverse 
temperature f3 if for any a, b £ A there exists a bounded holomorphic function 
F a j, continuous on < Im(z) < f3, such that, for all isl, 

F a , b (t) = (p(aat(b)) and F a>b (t + i/3) = <p(a t (b )a) (9) 

Equivalently, the KMS condition © is expressed as the relation 

<p(<r t (a) b) = <p(b<Tt+ip(a))- (10) 

Consider now a Cuntz-Krieger C*-algebra Oa- Any element in the *- 
algebra generated algebraically by the Si can be written as a linear com- 
bination of monomials of the form S^S*, for multi-indices /i = (ix>---iui) 
and v — (ji, . . . , The subalgebra Ta is the AF-algebra generated by the 
elements of the form S^S*, for = \i>\. It is filtered by finite dimensional 
subalgebras !Fk, for k > 0, generated by the matrix units Ei v = S^PiS*, with 
|^| = \u\ = k, where Pi = SiS* are the range projections of the isometries Si. 
The commutative algebra of functions on the Cantor set Aa of the subshift 
of finite type (Aa, T) associated to the Cuntz-Krieger algebra, is a maximal 
abelian subalgebra of Ta, identified with the elements of the form S^S* (cf. 
0)- 

In our case, one can consider a time evolution on O a n of the type consid- 
ered in [7], with a™' h (Sk) = exp(it(u—h)) Sk, for u = log rjp = P((3), the topo- 
logical pressure, i.e. the top eigenvalue of iJSJ, and h(x) = — (3/2 log|T'(x)|. 
Thus, for all t £ K, the function exp(-ith) acts on the elements of Oa n by 
multiplication by an element in C(Em x P). 

Then a KMSi state ip for this time evolution satisfies the relation (cf. 
Lemma 7.3) 

fc 

for ae Oa n - For all a = / e C(E N x P), we have £ fc S* k e h f S k = C h (f), 
where the Ruelle transfer operator 
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(9,r)eT-i(i, 8 ) 

is in fact the Perron-Frobenius operator JSJ). In particular, the KMS condition 
implies that the state ip restricts to C{En x P) to a probability measure \i 
satisfying L* h \i — e u fi. For u — lognp, the existence and uniqueness of such 
measure can be derived from the properties of the operator JHJl , along the lines 
of [12 0. 

Theorem 5.2 For (3 < 2\ogr(Aff)/ \og(N + 1), there exists a unique KMS\ 
state for the time evolution erf ^ 2 log ' T ' on the algebra Oa n - This re- 
stricts to the subalgebra C{En xP) to/H J f dfi with the probability measure 
satisfying C^n = r\ppi, for the Perron-Frobenius operator (jSJ. 

Proof. Since by Lemma 15.11 the matrix An is irreducible and aperiodic, by 
Proposition 7.6 of [7], there is a surjective map of the set of KMS states 
to the set of probability measures satisfying = e u \i. Uniqueness fol- 

lows from Lemma 7.5 and Theorem 7.8 of |7j, by showing that the estimate 
varo(h) < t{An) holds, where r(AN) is the spectral radius of the matrix An 
and varo (h) = max h — min h on En X P. This provides the range of values of j3 
specified above, since on En x P we have varo(— (3/2 log \T'\) — \og(N + \)P/ 2 . 
The KMS state ip is obtained as in [7], by defining inductively compatible 
states 

<p k {a) =e- u Y,Vk-i{S*e h / 2 ae h ' 2 S j ) for a e T k . 

3 
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